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§ . Latest results on structure functions, as available at the Lepton-Photon Symposium 

1999, are presented. This report focusses on three experimental areas: new structure 
function measurements, in particular from HERA at low x and high Q 2 ; results on 
light and heavy flavour densities; determinations of the gluon distribution and of a s . 
As the talk was delivered at a historic moment and place, a few remarks were added 
recalling the exciting past and looking into the promising future of deep inelastic 
scattering (DIS). 

> : 

5} : 1 Introduction 

o : 

About three decades ago, highly inelastic electron-proton scattering was observed 
q . by a SLAC-MIT Collaboration [|IJ which measured the proton structure function 

vW 2 {Q 2 \v) to be independent of the four-momentum transfer squared Q 2 at fixed 
Bjorken x = Q 2 /2M p u. Here v = E - E' is the energy transferred by the virtual 
V ' photon. It is related to the inelasticity y through v = sy/2M p , with proton mass M p 

and the energy squared in the centre of mass system s = 2M p E. With the SLAC 
linear accelerator the incoming electron energy E had been successfully increased 



by a factor of twenty as compared to previous form factor experiments 0. Thus 
Q 2 = AEE' sin 2 (6 l /2) could be enlarged and measured using the scattered electron 
energy E' and its polar angle 6. Partonic proton substructure was established at 
l/y/tp ~ 10~ 16 m which allowed the scaling behaviour fij of vW 2 {Q 2 ,v) — > F 2 (x) 
to be interpreted. In the quark-parton model (QPM) the structure function F 2 is 
given by the momentum distributions of valence and sea quarks, q = q v + q s , and of 
antiquarks q weighted by the square Q 2 of the electric charge, F 2 (x, Q 2 ) = x J2 q Qq(q+ 
q). Neutrino experiments found o v ~ 3 a D demonstrating that partons could be 
identified with quarks having gauge couplings like leptons and that at large x the sea 
is small. Scaling violations were hidden in the first DIS data taken at x ~ 0.2, as if we 
needed help to understand the basics of inelastic scattering. They were found in fiN 
scattering in an extended x, Q 2 range. The logarithmic Q 2 dependence of F 2 (x, Q 2 ), 
established in subsequent neutrino and muon-nucleon scattering experiments, was 
attributed to quark-gluon interactions in Quantum Chromodynamics |7J. 



1 



I 

H 



o.,s 

0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
O.J 

O.J 
0.2 











m BCDMS 1981 




— CTEQ4M 40 GeV 2 




f --- CTEQ4M 180 GeV 2 






> 


- 40<Cf<180 GeV 2 











O.J 0.2 0.3 0.4 0.5 0.6 0.7 0.1 



Figure 1: Measurement of the jZ interference structure function xG% in y u ± Carbon 
scattering by the BCDMS Collaboration compared with a recent parametrization of 
the valence quark distributions by the CTEQ group. 



With the discovery of neutral currents pi DIS neutrino experiments made a major 
contribution to the theory of weak interactions. In 1979 another ep scattering exper- 
iment was performed at SLAC |J which determined in a highly sensitive polarization 
asymmetry measurement at Q 2 ~ 1.5 GeV 2 the right-handed weak isospin charge of 
the electron to be zero. This experiment selected thus the Glashow Weinberg Salam 
model as the standard electroweak theory and opened the possibility to investigate 
proton structure at high Q 2 via Z boson exchange. The nucleon structure function 
F 2 was generalized, still in a V — A theory [[HJ, to three functions 

(F 2 , G 2 , H 2 ) = x J2(Ql 2 Q<,v q , v\ + a 2 q )(q + q) (1) 
q 

arising from photon exchange {F 2 ), 'yZ interference (G 2 ) and Z exchange (H 2 ), where 
v q (a q ) are the vector (axial vector) quark couplings In charged lepton-nucleon 
neutral current (NC) scattering two further structure functions appear which are 
analogous to xF% in neutrino scattering 

(xG 3 ,xH 3 ) = 2x^2(Q q a q ,v q a q )(q - q). (2) 
<? 

A DIS muon experiment with simultaneous beam charge and polarity reversal resulted 
in the first determination of the 7Z interference structure function xG% at Q 2 ~ 
60 GeV 2 by the BCDMS Collaboration at CERN, Fig. |. Electroweak interference 
occurs at the level of k ~ 10~ 4 Q 2 / GeV 2 as defined by the ratio of the weak and 
the electromagnetic coupling constants. Since the axial vector couplings could be 
considered to be known this was an interesting measurement of the valence quark 
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distribution sum u v + d v which confirmed the sign of the quark charge combination 
Qu ~ Qd to be positive. 

With the HERA energy of s = AE e E p ~ 10 5 GeV 2 the kinematic range of DIS 
experiments could be greatly extended towards high Q 2 since s was enlarged by a 
factor of about 2E P / GeV compared to fixed target scattering. The first measurements 
of F 2 by the HI |T| and the ZEUS @ Collaborations, using data taken in 1992, 
reached x ~ 0.0005 at Q 2 ~ 20 GeV 2 . They discovered a steep rise of F 2 (x,Q 2 ) 
towards low x at fixed Q 2 : below x ~ 0.01 a decrease by one order of magnitude 
translates into an increase of F 2 by about a factor of two, Fig. ^. Although a "Possible 



3.5 



2.5 



1.5 



0.5 



Q 2 =1 5 GeV 2 



®H1 92 
O ZEUS 92 




10 



10 



10 



1 

X 



3.5 



2.5 



1.5 



0.5 




10 



10 



10 



1 
X 



Figure 2: Measurements of the proton structure function F2(x,Q 2 ) by the muon- 
proton scattering experiments BCDMS and NMC and their extension towards low x 
by the first available HERA data on F 2 . The curve represents the anticipation of this 
rising behaviour by Gliick, Reya and Vogt fT7|j. 



Non-Regge Behaviour of Electroproduction Structure Functions" 15 at low x had 



been considered and the concept [16| and modified phenomenology |1_7] of 'dynamical 
partons' had been worked out, this rise came as some surprise since the DGLAP 
evolution equations do not a priori fix the x behaviour. This rise is now basically 
understood as being due to the dominance of gluons which leads to the description of 
the scaling violations as dF 2 /d\n Q 2 oc a s -xg for Q 2 3> M 2 and low x. Its quantitative 
description in NLO QCD and the search for new dynamics [18| connected with large 
logarithms of 1/x requires highest possible precision, i.e. improved instrumentation 
and higher luminosity than was available when the first observation was made. 

While much attention has been paid to the inclusive and charm structure function 
measurements at HERA, remarkable progress was also achieved in the investigation 
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of up, down, strange and charm quark distributions with neutrino and Drell-Yan 
experiments at the Tevatron. 

This paper describes a talk on structure functions in deep-inelastic scattering 
delivered in 1999. Such a report is to some extent personal and cannot possibly cover 
this expanding field of particle physics in any exhaustive fashion. It thus may be seen 
together with further articles, e.g. |19|, |20fl , and with the conference on deep inelastic 
scattering and QCD held at Zeuthen in April 1999 pi] . It demonstrates remarkable 
progress in DIS since the previous Symposium on Lepton-Photon Interactions ||22|| . 
This talk focussed on recent measurements of structure functions (Section 2), of quark 
distributions including charm (Section 3) and determinations of the gluon distribution 
and of a s (Section 4). The field of deep inelastic lepton-nucleon scattering has an 
exciting future as will be described briefly in Section 5. 



2 Recent Measurements of Structure Functions 

Since the first SLAC experiment, fixed target muon and neutrino-nucleon scattering 
experiments and subsequently the HERA collider experiments HI and ZEUS ex- 
tended the explored kinematic region of DIS by several orders of magnitude, Fig. |3|. 
At smallest x partons carry only a vanishing fraction of the proton momentum. Hence 
the kinematics resembles the fixed target experiments where both the electron and 
hadrons are scattered into the lepton beam direction (unfortunately termed 'back- 
ward' at HERA). For high Q 2 > sxE e /(E e + xE p ), i.e. Q 2 > 2, 800 GeV 2 for x > 0.5, 
the electron is scattered through angles 9 > 90° with respect to the electron beam 
direction, similar to Rutherford backscattering. The kinematic range of the HERA 
collider experiments is confined to about y > 0.001. For lower y hadrons escape in 
the forward (proton beam) direction. At very small y the inclusive kinematics can- 
not be reliably reconstructed using the scattered electron variables alone since the x 
resolution varies like 1/y. 

Until 1997 HERA ran with positrons scattered off protons of 820 GeV energy and 
about 40 pb _1 of luminosity became available for each collider experiment. From 1998 
till May 1999 data samples of about 15 pb _1 were collected in collisions of electrons 
with 920 GeV protons. The energy is tuned to about 27.5 GeV to optimize the 
polarization for the fixed target experiment HERMES. Longitudinal lepton beam 
polarization is foreseen to be used in colliding beam mode from 2001 onwards. 

2.1 Transition to Photoproduction and Low Q 2 < M 2 

The structure function F 2 which dominates the DIS cross section behaves like x~ x ^ x '^ 
and vanishes due to gauge invariance with Q 2 — > like 0(Q 2 ). The total virtual 
photon-proton scattering cross section is related to F 2 as o2 t — ^ 2 a ■ F 2 /Q 2 . Mea- 
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surements of F 2 at low Q 2 investigate the dynamics of the transition from the deep 
inelastic to the photoproduction regime [p3fl . In Regge theory the structure func- 



tion F 2 results from a superposition of exchanged Regge poles with intercepts a^, 
F2 = E j3i(Q 2 )W 2ai ~ 2 , where W 2 ~ Q 2 /x > Q 2 for low x, W being the invari- 
ant mass of the 7*p system. A recent fit to F 2 data (DL98) is rather successful 
using three trajectories, i.e. ot\ = 1.08 for the soft pomeron, a 2 = 0.55 for a and 
/ exchange and 03 = 1.4 for the so-called hard pomeron p3}|. For Q 2 — > the 



exponent A is approximately given by the dominant pomeron Regge trajectory, i.e. 
A ~ O! - 1 ~ 0.1. The recent ZEUS data PSI], obtained with a backward calorime- 



ter and tracker positioned close to the beam pipe, are rather well described by this 
model, Fig. [|. Phenomenological models using a combination of Generalized Vector 
Meson Dominance 0] and perturbative QCD |27| describe this transition also well. 
Extrapolations of F 2 (x, W 2 ) to Q 2 ~ come out to be somewhat higher than the 
direct measurements of <r7 oi p []25| with tagged electrons. The F 2 based a] t data are 
still at some Q 2 distance from the real photoproduction measurements which have 
uncertainties of about 10% due to beam optics and the imperfect simulation of the 
complete final state. Further extension of the range of the inclusive F 2 measurements 
at HERA towards lowest Q 2 values is thus desirable. This could be achieved in a 
rather short run of HERA at minimum possible electron beam energy since Q 2 is 
proportional to E 2 for all except the high y values. 
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Figure 4: Measurements of the proton structure function F 2 (x, Q 2 ) expressed as 
a tot(Q 2 i W 2 ) from recent low Q 2 data by ZEUS (closed points) and rebinned F 2 data 
at larger Q 2 by HI and ZEUS (open squares). The two triangles at W = 172 GeV 
and W = 189 GeV near Q 2 rs denote the direct photoproduction cross section 
measurements of HI and ZEUS. 



New data on parton-hadron duality p8| became available this year from an exper- 
iment at Jefferson Laboratory measuring electron-proton and deuteron scattering 
in the resonance region W 1 GeV. The superposition of cross sections, determined 
at different Q 2 between 0.2 and 3.3 GeV 2 , leads to an averaged behaviour of F 2 which 
is valence like even at low x, or mass corrected £ ]3(J, which supports the assumption 
made in the GRV analysis JT7| for the initial x distributions at very small Q 2 . In 
this experiment, which in the future will measure the ratio R = (Jl/&t, one estimates 
power corrections ('higher twists') to be small and derives the magnetic elastic proton 
form factor G P M from inelastic data. 



2.2 Neutrino Experiments 

The final measurement of vFe and vFe scattering cross sections by the CCFR Col- 



laboration 31 



is in good agreement with previous data obtained by the CDHSW 
Collaboration and more accurate. The high statistics CCFR data has been used for a 
number of investigations regarding all structure functions involved (Sections 
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and also for tests of QCD (Section |4.3| ). Recently data were released for extremely 
large x > 0.75 pointing to cumulative effects beyond Fermi motion in the nucleus |]32| 



which were studied previously by the BCDMS Collaboration [33 



Data were obtained by the IHEP-JINR neutrino experiment in the wide band 



neutrino beam at the Serpukhov U70 accelerator [33]. Based on about 750 v and 
6000 v events for W 2 > 1.7 GeV 2 and Q 2 ~ 2 GeV 2 , the structure functions F 2 and 
xF 3 were disentangled and a s {M 2 ) = 0.123l^ was determined in NLO QCD. 

2.3 Precision Measurement at Low x and Medium Q 2 

The HI Collaboration released for this conference the so far most precise measurement 
of the DIS cross section at HERA. In its reduced form it can be written as 

Q 4 x d 2 a y 2 

Or = F 2 - — ■ F L , (3) 



2vra 2 y + dQ 2 dx Y + 

i.e. o> ~ F 2 apart from high y where a r — > F 2 — Fl oc o>t- Here Fi denotes the 
longitudinal structure function which is related to the ratio R = F L /(F 2 — F L ) and 
Y + — 1 + (1 — y) 2 . The HI data, taken in 1996 and 1997, have statistical errors of 
typically 1% and systematic errors of 2-3%, apart from edges of the acceptance region. 
In order to reach this precision HERA has been anually increasing the luminosity. The 
HI experiment was subject to a major upgrade of its backward apparatus replacing a 
Pb-Scintillator calorimeter by a Pb-fibre calorimeter of higher granularity, an MWPC 
by a planar drift chamber and adding a high resolution Silicon strip detector telescope 
for electron track identification and kinematic reconstruction. This upgrade permitted 
the measurement to be extended to high y < 0.89 in order to access Fl (Section ^4|) 
and to low y > 0.003 in order to reach the x range covered by DIS fixed target 
experiments. Comparing the data shown in Fig. |5| with the initial HERA data, 
Fig. 0, one recognizes the impressive progress made in a few years. The data are well 



described by NLO QCD as discussed in Section (L2[ Consistent results on preliminary 
F 2 data were previously obtained by the HI and ZEUS Collaborations |J5] . 

The HI data help resolving a long standing controversy between NMC and E665 
HP data and the CCFR uN data on the structure function F 2 . As shown in Fig. |5] the 
HI data overlap and extrapolate well to the fip data. The CCFR F 2 determination 
which is being redone |3(| was recently criticized regarding the treatment of charm 



and shadowing P7| . Since F 2 and xF 3 add up to the measured cross section, an F 2 
reanalysis may affect also the value of a s derived from xF 3 . The CCFR cross section 
measurement improved in a consistent way the CDHSW cross section data. Those 



seem not to be in contradiction with muon data [pq| . 

Precision measurements at HERA are essential for calculating the expected rates 
at LHC energies and also permit to estimate the neutrino scattering cross sections in 
active galactic nuclei or gamma ray bursts at ultra high energies, up to E u ~ 10 12 GeV. 
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Figure 5: Measurements of the DIS cross section by HI compared with NMC and 
BCDMS up data. The solid curve denotes the fitted cross section in NL0 QCD using 
HI and NMC data for Q 2 > 3.5 GeV 2 . The dashed curve is the structure function F 2 
obtained in the QCD fit which at low x departs from a r . 



Recently very high energy rates were calculated using the DGLAP equations |9[ 



the GRV approach in DGLAP QCD g| and a combination of DGLAP and BFKL 



dynamics |4l]] which agree remarkably well. 



2.4 Longitudinal Structure Function Fl 

In the naive QPM the longitudinal structure function Fl is zero since partons have 
spin 1/2. In QCD it acquires a possibly large value due to gluon emission and repre- 
sents together with F 2 a strong constraint to the theory in NLO. 

The sum of v and v nucleon scattering cross sections is proportional to 2xF\{\ + 
eR)— Y_-AxF 3 /2Y + and thus is sensitive to R where e = 2(1— y)/Y + is the polarization 
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Figure 6: F 2 structure function data of CCFR, NMC, E665 and HI. The CCFR data 
were corrected for nuclear effects and for the difference of F 2 in charged lepton and 



neutrino scattering [36]. The CCFR data are shown with statistical errors only. 



of the W boson exchanged and Y± = 1±(1— y) 2 . The CCFR Collaboration has studied 
the Q 2 dependence of R for 0.015 < x < 0.5 and Q 2 < 5 GeV 2 using phenomenological 
descriptions for the strange and charm quark distribution difference determined by 
AxF 3 ~ Ax(s — c). The ratio R tends to be large, R > 0.5, at small Q 2 ~ 1 — 2 GeV 2 
and x < 0.1. For Q 2 > 10 GeV 2 the function AxF 3 = xF^ — xF^ was extracted which 



is of interest for the treatment of massive charm j42 . 

Using unpolarized targets the HERMES Collaboration measured the ratio of ni- 
trogen to deuterium electroproduction cross sections to be astonishingly small at low 
Q 2 [fi3"| . This effect has been attributed to a very large ratio Rn/Rd > 5 in the region 
0.01 < x < 0.06 and 0.5 < Q 2 < 1.5 GeV 2 with as yet unexplained origin. 

The measurements of the longitudinal structure function in ep and up scattering 
are summarized in Fig. [?[ The HI data were obtained using assumptions for the 
behaviour of F 2 in QCD (for Q 2 > 10 GeV 2 ) and, independently of QCD, for the 
derivative dF2/dhxy (for Q 2 < 10 GeV 2 ) in the high y region [34] where the cross 
section approaches F 2 — F^. Contrary to fixed target experiments such assumptions 
are possible since HERA covers more than two orders of magnitude in y where F 2 
can be fixed independently of F^. The overall behaviour of F^ as a function of x is 
well described by a QCD fit in NLO using F 2 data only, i.e. by deriving the gluon 
(and parton) distributions from scaling violations and then calculating Fl (Fig. [7]). 
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Figure 7: Measurements of the longitudinal structure function Fi(x,Q 2 ) in ep and 
HP scattering. At low x Fl is large because of dominant contributions of a large 
gluon momentum density. The four bins comprise data between Q 2 = 2, 5, 15, 25 and 
50 GeV 2 , respectively. The curves represent the HI calculation of F L . 



The behaviour of R observed at low Q 2 ~ 1 GeV 2 and the so far limited accuracy 
of the HI Fl data, obtained with 6.8 pb -1 , represent a challenge for forthcoming 
experiments and their theoretical interpretation. This comprises the hypothesis of 
particularly large higher twist effects and large higher order corrections which at low 
x and Q 2 may become even negative in NLO due to a large negative contribution of 
the gluonic coefficient function [ 45|] . 



2.5 Weak Neutral Currents at HERA 

At high Q 2 ~ Myyz photon, Z-boson and PU-boson exchange are of comparable 
strength. Thus electroweak interactions can be used to probe proton structure in 
neutral (NC) and charged current (CC) scattering at HERA in the same experiments. 
This is demonstrated with the Q 2 distributions in electron and positron proton NC 
and CC scattering, Fig. 
by ZEUS (e+ NC @, e 



j, measured by HI 
CC 
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The double-differential NC cross section, neglecting the three longitudinal struc- 
ture functions, is given by two generalized structure functions F 2 and xF 3 
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Figure 8: Measurements of the Q 2 dependence of the positron and the electron proton 
neutral and charged current scattering cross sections at HERA, using data taken in 
1994-1997 (e+ E p = 820 GeV) and in 1998-1999 (e~ E p = 920 GeV). Electromagnetic 
and weak interaction cross sections become of similar strength for Q 2 ~ M§, Myy- 



These depend on the quark couplings and distributions but, contrary to hadronic 
tensor definitions of structure functions fLOf , they depend also on the weak electron 
couplings v, a to the Z boson, on the longitudinal electron beam polarization (A) 
and on the propagators via k = Q 2 /[4 sin 2 6 W cos 2 6w{Q 2 + M%)] where 0w is the 



electroweak mixing angle. They comprise five genuine structure functions [] 1 



F 2 + k(-v =f \a)G 2 + k 2 (v 2 + a 2 ± 2Xav)H 2 



xFg = k(Xv ± a)xGz + k z 



-X(v 2 + a 2 ) T2av)xH : 



.3- 



(5) 
(6) 



defined in Section 1, Eqs. [I] and 0. The xF 3 term (oc Y_) contributes sizeably only 
at large y and high Q 2 . The high Q 2 NC cross sections measured currently at HERA 
for A = are approximately given by 



Y,-Fo± naY_ ■ xG^ 



(7) 



This causes a positive charge asymmetry between electron and positron scattering 
which is proportional to aa q , i.e. parity conserving, and which is determined by the 
function xG$ measured previously by BCDMS at lower Q 2 for an isoscalar target, see 
Fig □ 

The HI Collaboration has performed measurements of double differential NC 
scattering cross sections j|7| using 35.6 pb _1 of e + data |^6[ taken in 1994-97 at 
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E p = 820 GeV and 15 pb" 1 of e~ data (47| taken in 1998-99 at E p = 920 GeV. A 
comparison of the cross section measurements with electrons and positrons is illus- 
trated in Fig. [| which agrees with expectation based on the jZ interference in NC 
scattering. 
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Figure 9: Measurements of the double differential NC e + and e~ proton scattering 
cross sections by HI and NMC and BCDMS data. The high Q 2 HI data are well 
described by a QCD fit, which may even be restricted to Q 2 < 150 GeV 2 f47|, repre- 
senting a remarkable confirmation of the DGLAP evolution in Q 2 . A fit to solely HI 
and NMC data passes through the HI points but undershoots the BCDMS data at 
largest x significantly. 
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3 Light and Charm Quark Distributions 



3.1 Charged Currents and Up and Down Quarks 

New information on the up and down quark distributions became available from 
improved measurements of the charged current cross section at HERA by HI and 
ZEUS. The double-differential CC scattering cross section is given as 

- - • ( M " W ) 2 • • \Y W ± =F Y xW±] (8) 

dxdy~2n \Q 2 + M 2 W ) * 2 i Y +W 3 TY-xW s \ (8) 

where G is the Fermi constant and My/ the mass of the W boson. The CC cross section 
contains two structure functions for a given lepton beam charge and is proportional to 
s. The HERA energy is equivalent to 53.9 TeV neutrino beam energy in a neutrino- 
nucleon fixed target experiment. The energy dependence is damped for Q 2 > Myy- 
In the QPM the CC structure functions are combinations of up and down quark 
distribution sums, i.e. W 2 + = D + U, = U + D, xW£ = D-U and xW 3 ~ = U-D 
with U — x(u + c) and D = x(d + s). At large x > 0.3 the valence quark distributions 
u v and d v dominate the interaction cross sections, i.e. 

a(e + p^uX) oc U + (1 - y) 2 D -> (1 - ufxd v (9) 
a{e~p ^vX) oc U + (1 - y) 2 D -> xu v (10) 

for x — > 1. A complete set of double differential e^p CC cross section data was 
presented by HI using 36 pb _1 of positron-proton data (1994-1997) |46| and 15 pb -1 
of electron data (1998-1999) f47j. The U dominated e~p cross section was found 
to be about 5 times larger than the e + p cross section at Q 2 ~ 10, 000 GeV 2 . The 
e + p CC data of HI are consistent with the published measurement of the ZEUS 
Collaboration ^| based on 47.7 pb -1 , Fig. |T0[ The NC and CC measurements at 
high Q 2 are of particular interest for the determination of the d/u ratio at high x 
because their interpretation is free of nuclear corrections. Yet, an order of magnitude 
increase in luminosity is still required to access the high x region which represents 
one of the goals of the HERA luminosity upgrade programme. 

Deuterium binding corrections were recently reconsidered, and d v was adjusted to 



be larger than previously assumed [51], the ratio d v /u v for x — > 1 tending to 0.2. An 
enlarged d quark distribution fits to the charge asymmetry data in pp collisons. 
Violation of u and d quark symmetry in protons and neutrons, however, which was 
suggested to explain the difference between the CCFR and NMC F 2 data [p2]| , leads 



to too large a W asymmetry [53 
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Figure 10: ZEUS measurement of the CC e + p scattering cross section compared 
with U and D quark distributions, see Eq. [TU], and different QCD fits. 



3.2 Sea Quarks 

Interesting data become available on the flavour asymmetry in the nucleon sea. From 
a high statistics measurement of Drell-Yan muon pair production in pp and pd col- 
lisions at the Tevatron, the E866/NuSea Collaboration obtained for Jq(u — d)dx a 
value of -0.118 ± 0.011 at (Q 2 ) = 54 GeV 2 This confirms and also significantly 
improves the previous NMC result of —0.15 ± 0.04 which was derived from a mea- 
surement of the Gottfried sum rule fo[{F£ - F£)/x]dx = 1/3 + 2/3 • - d)dx. The 
measured ratio d/« as a function of x is shown in Fig. [IT]. The data have considerable 
impact on global parametrizations of parton distributions. A consistent result, albeit 
of less statistical accuracy, was obtained by the HERMES Collaboration with 
a measurement of semi-inclusive 7r ± production in unpolarized ep and ed scattering 
at lower (Q 2 ) = 2.3 GeV 2 . A violation of flavour symmetry is not predicted in per- 
turbative QCD which points to non-perturbative effects such as Pauli blocking and 
pion clouds. In the latter model the nucleon is expanded in a Fock state of mesons 
and baryons. Phenomenologically one finds more 7r + than 7r~ in the nucleon with a 



momentum distribution peaking at x n ~ 0.2 [56 
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The NuTeV Collaboration []57|] determined the strange quark distribution to be 
about 1/2 of the averaged nucleon sea, i.e. s = [0.42 ± 0.07(syst) ± 0.06(stat)] ■ {u — 
d)/2, in agreement with previous analyses of dimuon production in neutrino-nucleon 
scattering experiments. 

Indications for a difference of the strange and anti-strange quark distributions at 
large x ~ 0.6 were obtained in a recent reanalysis and global fit of DIS and Drell-Yan 
|38[ |. Sensitivity to (s — s) in this analysis comes from the CDHS data measuring 
a u oc x(s — s) + Y^x{u v + d v ) at high x. Such a strange asymmetry is possible 



data 



a 



in models considering states as K + A to be intrinsic to the nucleon where K + yields 
s oc (1 — x) and A yields s oc (1 — x) 3 
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3.3 Charm 



Charm, as was already noticed by Witten in 1976, may "subject non-Abelian theories 
to a rigorous experimental test by measuring the charmed quark contribution to 
structure functions" |5]| . Since then the charm and beauty treatment in perturbative 
QCD has been worked out to higher orders Variable flavour schemes are being 



studied |61[ to correctly handle the heavy flavour contributions near and beyond 
threshold in analyses of parton distributions, of the gluon distribution and of a s . A 
new measurement of the charm structure function F 2 CC was performed by the ZEUS 
Collaboration [62] using the AM tagging technique for D* —>■ K2tt and K4n, Fig. |I~2| . 
The relative contribution of charm is large, reaching 30% at low x < 0.001 for Q 2 ~ 
100 GeV 2 . This large fraction is due to photon-gluon fusion as the dominant process 
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Figure 12: Measurement of -F| c (x, Q 2 ) and of the ratio / F 2 in ep scattering at 
HERA by the ZEUS Collaboration with 37 pb -1 of data. The dashed error bands 
denote the uncertainty of the QCD fit which is dominated by the charm quark mass 
range chosen to be 1.2 to 1.6 GeV. 



for charm production. Further experimental progress at HERA towards high precision 
will be achieved with new or upgraded Silicon vertex detectors, higher luminosity, 
inclusion of further final states and dedicated track triggers. 

4 Gluon Distribution and Coupling Constant ct s 
4.1 Scaling Violations at Low x 

Scaling violations in the DIS Q 2 region down to low x ~ 0.00005 can be successfully 
described in the DGLAP formalism. This is again demonstrated with the new precise 
cross section measurement of HI, Fig. |5|. Conventional QCD fits use parametrizations 
of parton distributions at a starting scale Q 2 Q and evolve them in Q 2 to highest Q 2 > 
M§ values up to order a 2 . However, the splitting functions have expansions which 
contain also powers of ln(l/a;). These are large at low x, such that a s ln(l/a;) ~ 1, 
and yet do not seem necessary to phenomenologically describe the observed structure 
function behaviour. Calculations are performed in order to account for these ln(l/x) 
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terms [20| and to cure perhaps the instability of the BFKL equation in NLO |63" |. 
Indications were reported for the presence of ln(l/x) terms in inclusive DIS data ]64 . 
Experimentally even higher precision is both required and possible for the structure 
function measurements, including Fl, which may lead to crucial tests of QCD at low 
x. Due to unitarity constraints one expects to find saturation of the rising behaviour 
of F 2 which, however, seems to be beyond the low x range accessible by HERA in the 
DIS region. 

Scaling violations are conveniently studied using the InQ 2 derivative of F 2 . In 



Fig. [13] the structure function F 2 from HI is shown as a function of Q 2 for x < 
0.01. The InQ 2 dependence is non-linear and can be well described by a quadratic 
expression P 2 — a + blnQ 2 + c(lnQ 2 ) 2 (solid lines) which nearly coincides with the 
NLO QCD fit (dashed lines). The local derivatives dF 2 /dhiQ 2 determined from the 
new HI F 2 data are not constant in Q 2 and also depend on x. Approximately they 
can be described for each bin of x by b + 2 • clnQ 2 . Small deviations from this 
behaviour occur in NLO QCD. Using this expression the derivatives are determined 
at fixed Q 2 and displayed as functions of x in Fig. [TBI. There is no departure observed 



from a rising behaviour of the InQ 2 derivatives down to Q 2 = 3 GeV 2 . If such a 
plot is made as a one-dimensional distribution, using the derivatives calculated for 
each bin of x at the mean Q 2 of a given bin, then the derivative dF 2 /d InQ 2 flattens 
starting at Q 2 ~ 6 GeV 2 [pq| . In the region covered by the HI data this behaviour 
reflects the restriction of the kinematic range of the measurement. Some analyses of 
the ZEUS data extending to lower Q 2 ~ 1 GeV 2 introduce screening corrections in 
order to describe the behaviour of F 2 pB| , |67[| . Both F 2 and Fl in this region should 
be measured with still higher accuracy (see Section 2.1.) as these permit important 
information to be deduced on the dynamic interplay of gluon and sea distributions, on 
the effect of higher order and power corrections and on the shadowing phenomenon. 

4.2 Gluon Distributions 

In QCD the Q 2 evolution of F 2 is governed by the strong interaction coupling constant 
a s . The evolution relates the quark distributions to the gluon distribution xg. The HI 
Collaboration has performed a new NLO QCD fit to the HI and NMC inclusive cross- 
section data. It uses the DGLAP evolution equations for three light flavours with the 
charm and beauty contributions added according to the NLO calculation of the boson- 



gluon fusion process p8fl . The proton structure function F 2 is a superposition of two 
independent functions with different evolutions, i.e. F 2 = 5/18 • S + 1/6 • N, where 
the singlet function S — U + D is the sum of up and down quark distributions and 
the non-singlet function iV = U — D is their difference. In the new HI fit a different 
linear combination is introduced such that U = 2/3 ■ V + A and D = 1/3 ■ V + A. 
In a simplified parton model ansatz with u = d and s + s = (u + d)/2 one finds 
V = 3/4 ■ (3u v — 2d v ). This allows the quark counting rule to be applied which 
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Figure 13: New preliminary HI data show F 2 (x,Q 2 ) to be non-linear in InQ 2 at 
low x (left). The derivative dF 2 / d\n.Q 2 is a continuously falling function of x for 
Q 2 > 3 GeV 2 (right). 



constrains J Vdx = 3. This ansatz is used to fit the cross-section data, Fig. |5], for 3.5 < 
Q 2 < 3000 GeV 2 assuming a s [M\) = 0.118. It is written in the MS renormalization 
scheme and generalized to account for the measured difference u—d and the fraction of 
strange quarks, see Section ||[ The salient feature of this new analysis is that it applies 
to DIS proton data only but correctly determines the gluon momentum fraction to be 
about 0.45 at Q 2 = 10 GeV 2 . The gluon distribution resulting from this fit is shown 
in Fig. [TJ| (left). The inner error band defines the experimental uncertainty of a few 
per cent at low x using the treatment of correlated systematic errors of []69[] . The 
outer error band comprises uncertainties due to dependencies on the fit parameters 
(.QmimQoi a si m c) an d on the choice of parametrizations for the initial distributions. 
A remarkable feature of xg is the crossing point at x ~ 0.06 which is analogous to 
the Bjorken scaling behaviour of F 2 and reflects the conservation of the gluon and 
quark momenta. In Fig. [14] (right) the gluon distribution is seen to agree very well 
with xg unfolded from the charm structure function DIS and photoproduction data 
of HI which confirms hard scattering factorization. It has early been recognized 
that in photoproduction (Q 2 ~ 0) the charm mass provides a hard scale |fTT |. 

While xg at low x is well determined by the HERA structure function measure- 
ments, there are sizeable uncertainties of one order of magnitude at high x ~ 0.6 
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Figure 14: Determination of xg by HI using NMC and HI Ip data in NLO QCD 
(left). Comparison of xg from scaling violations with the unfolded gluon distribution 
from charm D* measurements by HI (right). 



The gluon distribution is accessed at high x by quark-gluon Compton scattering 
leading to direct photon emission [173]. In a recent experiment by the E706 Collab- 
oration WM the photon p T spectrum was found to exceed QCD expectation by a 



factor of about two which has been phenomenologically cured by a Gaussian trans- 
verse momentum smearing with kx of 1 GeV, larger than the intrinsic kx value of 
about 0.4 GeV |7^| . High E T jet data at large rapidities are sensitive also to xg at 
large x and lead to a rather high gluon distribution. Inclusion of different data sets 
yields remarkably differing results. Resolving the issue of xg at high x is essential 
for a reliable prediction of Higgs production in pp colliders. It is necessary since the 
high x exponent c g of xg oc (1 — x)° 9 is known to be correlated with a s . In this 
respect precision measurements of structure functions at high x are important. Since 
F 2 vanishes as (1 — x) 3 , any measurement error at large x is amplified like 1/(1 — x). 
The HERA collider experiments with their unique possibility to overconstrain the 
kinematics can be expected to lead to precision data also at high x J75[] when the 
luminosity is upgraded. 

Recently updates of the GRV parametrizations were presented f76|| . New sets of 
fits were made by the MRST |77[ and the CTEQ groups |78j. GRV98 uses DIS, n/p 
and Drell-Yan data assuming a s (M|)=0.114. MRST99 uses direct photon data of 
the WA70 experiment for different kx and varies the d/u ratio, a s and m c . CTEQ5 
does not use direct photon data but analyzes high Ex jet data instead. Sets are pro- 
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vided for different renormalization schemes and heavy quark treatments. As a conse- 
quence there exists a variety of parametrizations illustrating the still large flexibility 
of theoretical assumptions and pointing to possible experimental contradictions. An 
interesting attempt was made recently [|79| to quantify the experimental uncertainties 
of parton distributions resulting from global QCD fits to DIS data. 



4.3 Determinations of a. 



s 



New determinations of cn s (Mf) with structure function data were presented recently. 
Conventional analyses parametrize a set of input quark distributions and xg at certain 
input scale Q 2 using the DGLAP equations to NLO to calculate the theoretical ex- 
pectation. Minimization of a x 2 function determines a s and the roughly 10-15 parton 
distribution parameters. The treatment of systematic errors affects both the central 
value and the error size of a s (Ml). At low Q 2 power corrections to the logarithmic 
evolution may be sizeable and anticorrelate with a s . Since analyses differ in these 
assumptions and use different sets of data, one may not be surprised to still find some 
spread of the quoted values of a s (M|). Using the SLAC, BCDMS and NMC p and 
n structure function data and taking into account systematic error correlations and 
higher twists oc 1/Q 2 , a value of a s (M 2 ) = 0.1183 ±0.0021 (exp) ± 0.0013(t%) has 
been derived ||80|| . A similar analysis |8l] including the published HERA data and 
adding all errors in quadrature yields a s (M|)= 0.114 ±O.OO2(ea;p)ig;0O4(i/M/) which is 
closer to a previous determination of a s (M§) based on SLAC and BCDMS data [82 |. 



The quoted theoretical errors represent the uncertainties of the renormalization scale 
fi r , the former analysis compensating part of the \x r dependence with the higher twist 
contribution. 

The theoretical uncertainties are diminished in NNLO calculations. So far only 
partial results are available on the 3-loop splitting functions while the (3 function 
and the coefficient functions are known [[83 1. This gave rise to a revival of moment 



analyses. In [84| the xF 3 data of the CCFR Collaboration are reconstructed us- 
ing orthogonal Jacobi polynomials. Power corrections are considered and a value of 
a s (M|)= 0.118 ± 0.002 (stat) ± 0.005 (syst) ± 0.003 (thy) is obtained in NNLO 
corresponding to 0.120 in NLO. While this uses a pure non-singlet function, not cou- 
pled to the gluon distribution, a new analysis of SLAC, BCDMS, NMC, ZEUS and 
HI data using Bernstein polynomials of F 2 yields a s {M 2 z ) = 0.1163 ±0.0023 in NNLO 
with a single error supposed to comprise all experimental and theoretical uncertain- 
ties. This analysis |85| is extended to a Q 2 range of 2.5 to 230 GeV 2 and includes 



power corrections. Its NLO result is 0.1175, and moments of xg are determined. 

Although all these analyses represent quite remarkable theoretical and experimen- 
tal progress, one still has to be cautious. The systematic error treatments of these 
analyses differ. An important issue is the possible incompatibility of different data 
sets. For example, the combination of SLAC and BCDMS data yields an a s (M|) value 
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near to 0.114. Yet, this is known to result from a superposition of the BCDMS data 
favouring a value of about 0.110 with the SLAC data preferring a s ~ 0.120. Further- 
more, the moment analyses, while theoretically advanced to NNLO, shift the data 
weight to large x where the accuracy of the data is less impressive. Moreover, there 
is a dependence of the result on the minimum Q 2 considered |n[ which often leads to 
the introduction of power corrections with phenomenological x dependence. Finally 
the likely presence of ln(l/x) terms will affect the data interpretation. It is thus con- 
cluded that the great potential of DIS data to determine a s (Mf ) requires still much 
more work in order to determine a s at the one per cent level of accuracy. 

Interesting ideas are pursued to replace in the QCD analysis xg by the derivative 
dF 2 / d\n.Q 2 |81 , [86] and to develop the method of truncated moments |87[ in order to 



avoid the low x region in analyses of structure functions other than F 2 . The approach 
of double asymptotic scaling at low x of F 2 [15 represents a three parameter solution 
of the DGLAP equations and may lead to a particularly accurate determination of 
a s (M|) [^, j39| . This solution predicts a steady increase of xg towards low x which 
yet has to be damped at certain x and Q 2 since xg may not exceed the proton size 



iir 2 by too big an amount 



901 



5 On the Future of Deep Inelastic Scattering 



During the year 2000 the HERA luminosity will be upgraded |9T[ in order to provide 
an integrated luminosity of 150 pb _1 per year. Variations of proton and electron 
beam energies and the use of electron polarization in colliding mode will further 
enable the electroweak structure function measurements and enhance the discovery 
potential of the machine. The modifications of HERA are accompanied by major 
detector upgrades of the luminosity, forward tracking and Silicon vertex detectors of 
HI and ZEUS. 

The main injector neutrino oscillation detector at Fermilab (MINOS near ) will lead 
to precise, high statistics data (~ 4-10 7 events/year) on the six structure functions (F 2 , 
xF 3 and Fl for vFe and vFe scattering) which is necessary to disentangle the nucleon 
sea, i.e. to measure u + d, c and s [|9*2]] . Measurements of the nuclear dependence 
of neutrino DIS cross sections using additional targets will determine vA shadowing 



and perhaps help resolving the CCFR-NMC puzzle, Section |2.2| . Increase of neutrino 
energy by a factor of 10 would be possible in a 250 GeV muon storage ring providing 
extremely intense neutrino beams f)3|, |94| . 



The obvious next step in electron-proton DIS is a new ep machine p5 |. The 



proposed linear collider at DESY, TESLA, may provide collisions of electrons of 
up to about 500 GeV against HERA protons of nearly 1 TeV. A similar energy of 
y/s ~ 1.5 TeV can be obtained in ep collisions at LEP-LHC energies. These machines 
differ in technology, luminosity and kinematics. Yet one can envisage extending the 
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Figure 15: Kinematic region of a possible future ep collider using 500 GeV electrons 
from TESLA and 920 GeV protons from HERA. The line y = 0.05 represents the 
upper kinematic limit of ep collisions in HERA. Since the TERA machine is symmetric 
in energy, it provides full containment of highly energetic electrons and hadrons at 
large x and high Q 2 as indicated with the iso-# lines at 20° and 160°. 



low x acceptance by a factor of 20 and DIS data to Q 2 ~ 500, 000 GeV 2 and beyond, 
(Fig. [15D. Saturation and sub-structure will be searched for in this extended range. 

30 years after the pioneering SLAC ep experiments deep inelastic scattering still 
has an exciting future. 



6 Concluding Remarks 

HERA has opened the field of low x physics which is governed by gluon interactions 
and which is far from being fully understood. The gluon momentum density at low 
x is very large. This causes the structure function to rise at low x, it determines 
the longitudinal structure function to be large and the production cross section of 
heavy flavours to be sizeable. Increasing experimental precision leads to sensitive 
tests of QCD at higher orders perturbation theory. Most accurate simultaneous de- 
terminations are in reach of the gluon distribution and the strong interaction coupling 
constant with DIS data. Electroweak neutral and charged current structure functions 
provide new insights in the proton structure at high x. Measurements at Q 2 ~ M§ 
probe the proton nearly 100 times below the parton level reached three decades ago. 
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It is a spectacular result that no substructure of leptons or quarks has been observed 
so far. At the same time significant progress is made with various fixed target and 
pp experiments leading to deeper insight in the partonic structure of the proton. The 
gluon distribution at large x is small but remains to be determined. The next step is 
in reach for tests of the inner proton structure down to 2 • lCT 19 m. The outcome is 
unknown and deep inelastic physics therefore worth continuing effort. 



It is a pleasure to thank John Jaros and co-organizers for an excellent Sympo- 
sium. Many thanks are due to colleagues from the various DIS experiments for pro- 
viding information and guidance in understanding their results. I have to thank too 
many individuals to be named here, members of the HI Collaboration and its struc- 
ture function group, physicists and engineers of the Zeuthen Silicon detector group, 
colleagues in the DIS99 conference committees, many theoretical and experimental 
physicists for useful discussions and reading the manuscript and also several known 
physicists around the BCDMS Collaboration who introduced me to deep inelastic 
scattering years ago. Modern particle physics is a huge common effort of a large, 
mostly friendly community. Particular recognition is due to the youngest: I sincerely 
thank Vladimir Arkadov, Doris Eckstein, Alexander Glazov and in particular Rainer 
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